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Inclusive \ V u b\ measurements at BABAR 

A Sarti on behalf of the BABAR Collaboration 

Universita di Ferrara, Dipartimento di Fisica and INFN, 1-44100 Ferrara, Italy 

We present two inclusive measurements of Cabibbo-Kobayashi-Maskawa (CKM) matrix element \V u b\: one uses the lepton energy 
spectrum (£)) and the other uses the invariant mass of the hadronic system (m x ) to discriminate signal (S — > X u Cv) and background 
(B — » X c tv) events in B — > Xtv transitions. Both analyses are based on data samples collected by the BABAR detector at the PEP-II 
asymmetric-energy B Factory at SLAC. 



> : 
00 . 

in ■ 

: 

o . 
m ■ 

o ; 
^ : 
1 

^ ■ 



X 



The element \V„b\ of the CKM matrix [ 1| plays a cen- 
tral role in tests of the unitarity of this matrix: it's ex- 
traction, based on tree level decays, gives results that are 
independent of new physics contributions. We report the 
determination of \V u i,\ from two different measurements 
of the inclusive charmless semileptonic B branching frac- 
tion, &{B — > X u fv) 1 , using Ei spectrum (endpoint) and the 
mx spectrum on the recoil of fully reconstructed B mesons 
respectively. 

The selection of B — > X u (v events is hampered by the pres- 
ence of a large B — > X c iv background: E/ and mx spectra 
are used to discriminate the two different transitions. The 
endpoint analysis is sensitive to approximately 10% of the 
Ei spectrum while the acceptance for the mx approach is 
larger: 70% of the mx spectrum is selected by analysis 
cuts. The extrapolation of the measured rates to the full 
phase space introduces theoretical uncertainties 2 [ 2|. Re- 
sults also depend on the shape function (SF) modeling of b 
quark Fermi motion inside the B meson. 
Both measurements are based on data recorded by the 
BABAR detector [ 3 1 at the PEP-II asymmetric-energy e + e~ 
storage ring at SLAC. The endpoint analysis data sample 
consists of about 23 million BB pairs (21 ftr 1 ) collected at 
the Y(4S) resonance (ON -peak), with an additional sam- 
ple of 2.6 ftr 1 recorded about 40MeV below the 7(45) 
peak (OFF-peak), while the mx measurement uses a data 
sample of about 88 million BB pairs (82 ftr 1 ) ON-peak. 
Monte Carlo (MC) simulations of the BABAR detector 
based on GEANT 4 [ |4] are used to optimize selection cri- 
teria and to determine signal efficiencies and background 
shapes. To simulate B — > X c lv transitions three models 
are employed: B — » D*(v decay is modeled following a 
parametrization of form factors based on HQET [ 5|; for 
B — > Dlv decays and higher mass charm meson states 
B — > D^€v the ISGW2 model [|6| is used; nonresonant 
decays, B — > D^nlv, are modeled according to a prescrip- 
tion by Goity and Roberts [[7). In the endpoint analysis the 
MC simulation of B — > X u iv events is based on the ISGW2 



Charge-conjugate states are implied throughout this paper. 
2 The results presented are depending on the parton-duality to which no 
error is assigned. 



model: the hadrons X u are represented by single particles 
or resonances with masses up to 1.5GeV/c 2 and nonres- 
onant contributions are not included. In the m x analysis 
B — > X u lv transitions are simulated with an hybrid model 
which is a mixture of resonant and nonresonant compo- 
nents. The Fermi motion of the b quark inside the B meson 
is implemented in the nonresonant component using the SF 
parameterization described in [ 8 1, and the fragmentation is 
handled by Jetset7 . 4 [|5|. 

1 Endpoint analysis 

For this analysis, electron candidates are selected in the 
momentum range from 1.5 to 3.5 GeV/c in the T(4S) rest 
frame with a solid angle defined by the electromagnetic 
calorimeter acceptance. 

The inclusive electron spectrum for charmless semilep- 
tonic B decays, measured in the Y{4S) rest frame in the 
momentum range of 2.3-2.6 GeV/c, is used to extract 
S(B — > X u tv). To suppress low-multiplicity QED pro- 
cesses and continuum processes consisting of nonresonant 
e + e~ — > qq production (q = u, d, s, c) at least three charged 
tracks per event are required and a cut on the ratio of Fox- 
Wolfram moments H2IH0 < 0.4 f 1101 is applied. The 
missing momentum four- vector /?„„,, = p,- Pb„ cb — Px~ Pi, 
where all momenta are measured in the laboratory frame 
and pi refers to the four-momentum of the initial state of 
the colliding beams, can be used to select semileptonic 
events: \p m i SS \ is requested to be larger than 1 GeV/c, to 
point into the detector fiducial volume and the angle be- 
tween the electron candidate and the missing momentum 
is required to be greater than n/2. Candidate electrons are 
rejected if, when paired with an opposite-sign electron, the 
invariant mass of the pair is consistent with the J/tf/ mass 
(3.05 < M e+e - < 3.15 GeV/c 2 ). For the selection crite- 
ria described above, the detection efficiency for charmless 
semileptonic decays in the electron momentum interval of 
1 .5 - 2.6 GeV/c ranges from ~ 0.4 to ~ 0.25. 
The raw spectrum of the highest momentum electron af- 
ter the subtraction of continuum background (determined 
from OFF resonance data sample) is shown in Fig.[^. Also 
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Figure 1. Electron momentum spectrum in the T(45) rest frame, 
(a) ON-peak data after continuum subtraction (solid triangles), 
and MC predicted background from BB events (B -» X u ev) (open 
triangles), (b) ON peak data after subtraction of continuum and 
MC predicted B -» X„ev backgrounds (data points with statis- 
tical errors). For comparison, the histograms show the expected 
signal spectrum from B — > X„ev decays, (c) The differential rate 
S(B — > X u ev) as a function of the electron momentum. The data 
(statistical errors only) are compared to the prediction (solid line). 



shown are the MC predictions of the expected signal from 
B — > X u ev decays and background contributions from all 
other processes. The result of the subtraction of all back- 
grounds is shown in Fig. \l]p. For a given interval in the 
electron momentum, the inclusive partial branching ratio 
is calculated according to 



AS : 



NON ~ A^OFF - NB-»X„ev 

2eN R - R 



(l +<W)- 



(1) 



Here Non refers to the number of electrons detected ON- 
peak and jVoff refers to the fitted continuum background 
in a specified momentum interval, Nb-»x u hv is the back- 
ground from BB decays, e is the total efficiency for de- 
tecting a signal electron from B — > X u ev decays (including 
bremsstrahlung in detector material), and 6 mc t accounts for 
the distortion of the electron spectrum due to final-state ra- 
diation. As the overall normalization the total number of 
produced BB events is used, N BB = (22, 630+ 19+362)- 10 3 . 
The systematic error introduced by the efficiency estima- 
tion for the signal events is 5%. The uncertainty in the 
continuum background subtraction is 5%. The error com- 
ing from the BB background modeling translates to a rel- 
ative error of 3%. Variations of the colliding beam energy 
introduce a systematic error in the B — > X c ev background 
subtraction (5%). The total systematic error on the partial 
branching ratio measurement is ~ 9%. 
The fully corrected differential branching ratio as a func- 
tion of the electron momentum is shown in Fig.^. Inte- 
grating over the interval from 2.3 to 2.6 GeV/c, we get: 



AS(B -> X u ev) = (0.152 ± 0.014 + 0.014) • 10" 



(2) 



To determine the charmless semileptonic branching frac- 
tion S(B — > X u ev) from the partial branching fraction 



AS(Ap), the fraction f u {Ap) of the spectrum that falls into 
the momentum interval Ap is needed. The CLEO collabo- 
ration has recently used the measurement of the inclusive 
photon spectrum from b — > sy transitions [ II II to derive 
fu(Ap) for B — > X u ev transition. They quote a value of 
f u (Ap) = 0.074 + 0.014 + 0.009 for the interval Ap from 
2.3 to 2.6 GeV/c. Relying on the CLEO measurement, the 
result presented here translates into a total branching ratio 
of S(B -> X u ev) = (2.05 + 0.27 exp + 0.46 f J ■ 10~ 3 . We ex- 
tract \Vub\ from the measured inclusive charmless semilep- 
tonic branching fraction with the relation in [ 1151 and the 
average B lifetime of Tg = 1 .608 ± 0.012 ps [ 15 1 and find 

\V ub \ = (4.43 + 0.29 e , p ±0.25 o/J£ + 0.50 / „ + 0.35. sy )-10- 3 (3) 

Here the first error is the combined statistical and system- 
atic error, and the second refers to the uncertainty on the 
extraction of \V u b\ from relation in f 1151 . The third one is 
taken from the CLEO analysis and is related to the exper- 
imental determination of /„. The last error accounts for 
uncertainties related to assumption that b — > sy spectrum 
can be used for shape function modeling in B — > X u ev tran- 
sition. 



2 \V ub \ measurement using the recoil of fully 
reconstructed B mesons 

This analysis is based on BB events in which one of 
the B meson decays in a fully reconstructed hadronic fi- 
nal state (B mo ) and the other one is identified as decay- 
ing semileptonically by the presence of an electron or a 
muon. The full reconstruction of one of the two B mesons 
reduces the overall efficiency, but allows to reconstruct 
both the neutrino and the hadronic system (X), to deter- 
mine the flavour and to separate charged and neutral B 
mesons. In order to reduce systematic uncertainties due 
to efficiency determination we extract the branching ra- 
tio R u/sl = S(B -> XJv)/S(B -> XZv) after measuring the 
number of events with one identified lepton. 
To fully reconstruct a large sample of B mesons, hadronic B 
decays of the type B reco — > D w 7 are selected. Y represents 
a collection of hadrons with a total charge of + 1, composed 
of Witt* tizK* niK® ri4n°, where n\ + «2 < 6, «3 < 3, and 
«4 < 3. The kinematic consistency of a B reco candidate 
with a B meson decay is checked using two variables, the 

beam energy-substituted mass toes = Js/4 - /?J and the 

energy difference, AE — Eb - ~\fs/2. Here V* refers to 
the total energy in the T(45) center of mass frame, and ps 
and Eb denote the momentum and energy of the B reco can- 
didate in the same frame, respectively. In events with more 
than one reconstructed B decay, the decay mode with high- 
est purity is selected. 

Semileptonic B decays, B — > Xiv, recoiling against the 
B reco candidate are identified by an electron or muon candi- 
date with a minimum momentum (p*) greater than 1 GeV/c 
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Figure 2. a) Signal MC m x distributions with the requirement of Figure 3. The m x distribution in B — > Xlv decays, a) Data (dots) 

a lepton with p* > 1 GeV/c. Measured m x distribution before and and fit components, b) Background subtracted data and signal 

after all cuts. b)Fit to the otes distribution for the lepton sample MC. 
with p* > 1 GeV/c in the recoil of a B reco candidate 



in the B rest frame. Correlation between the charge of 
the prompt leptons and the flavor of the B reco is imposed 
(B° - B° mixing rate is used to extract the prompt lepton 
yield in case of neutral candidates). 

The hadron system X in the decay B — > Xlv is made 
of charged tracks and neutral energy depositions in the 
calorimeter that are not associated with the B reco candidate 
and not identified as a lepton. The mass of the hadronic 
system is determined by a kinematic fit that imposes four- 
momentum conservation, the equality of the masses of the 
two B mesons, and forces pL iss = 0. 
The selection of B — » X u lv decays is tightened by requir- 
ing exactly one charged lepton with p* > 1 GeV/c, charge 
conservation (Qx + Qe + Qb,„„ = 0), and missing mass 
consistent with zero (pf nlss < 0.5GeV 2 /c 4 ). These criteria 
improve the resolution in mx and suppress the dominant 
B — > X c £v decays, many of which contain additional neu- 
trinos or undetected Ki. We suppress the B° -> D* + Tv 
background with a partial reconstruction in which only the 
slow pion from the D* + — > Dn + S decay and the lepton are re- 
constructed. We veto events with charged or neutral kaons 
in the X system to reduce the background from B — > X c lv 
decays. The impact of the selection criteria on the mx dis- 
tribution is illustrated on MC in Fig.|5^. We determine R u / S i 
from N u , the observed number of b — > u events, and N s u 
the number of events with at least one charged lepton: 



hi ■ 



S(B X u lv) NJiV mx ) 



S(B -> Xlv) N sl 



^sl si 
sefm x > fc ; fc > 
x — ~ — ~~ . 



(4) 



Here s" , = (34.2 ± 0.6)% is the efficiency for selecting 
B — > X u £v decays with all analysis requirements, sJL = 
(73.3 + 0.9)% is the fraction of signal events with mx < 
1 .55 GeV/c 2 , sf Itf = 0.887 ± 0.008 corrects for the differ- 
ence in the efficiency due to the lepton momentum cut for 
B -> Xtv and B -> XJv decays, and efle'j = 1.00 ± 0.04 
accounts for a possible efficiency difference in the B reco re- 
construction in events with B — > Xlv and B — > X u lv de- 
cays. 



We derive N s i from a fit to the m E s distribution shown in 
Fig-Eb- The fit uses an empirical description r 1 1 31 of the 
combinatorial background from continuum and BB events, 
together with a narrow signal [ 14 1 peaked at the B meson 
mass. The residual background in N s i from misidentified 
leptons and semileptonic charm decays amounts to 6.8% 
and has been subtracted. We obtain N u from the mx dis- 
tribution with a x 1 fit to the sum of three contributions: 
signal, background A^ e from B — > X c lv, and a background 
of less than 1 % from other sources (misidentified leptons, 
secondary r and charm decays). 

In each bin of the mx distribution, the combinatorial B reco 
background is subtracted on the basis of a fit to the (Hes 
distribution. Figure^ shows the mx distribution with the 
results of the fit superimposed. The fit reproduces well the 
data having ax 2 /dof = 7.6/6. In the fit, the first bin is 
chosen to contain all events with mx less than 1.55 GeV/c 2 
while the other bins are chosen in order to separate the con- 
tribution from each resonant B — > X c (v mode. The mx cut, 
set at 1.55 GeV/c 2 , has been optimized minimizing the to- 
tal error. Figure |5p shows the mx distribution after back- 
ground subtraction with finer binning. Table ^ summa- 
rizes the results of fits with different requirements on mx, 
for electrons and muons, for neutral and charged B reco can- 
didates, and for different ranges of the B reco purity, V . The 
results are all consistent within the uncorrected statistical 
errors. 

We have performed extensive studies to determine system- 
atic uncertainties. We use events with charged and neutral 
kaons in the recoil of the B reco candidate as a control sam- 
ple to assess that the background from B — > X c £v events 
is properly described. The relative systematic error (A r ) 
due to the selection criteria related to the reconstruction 
of particles in the event is A r = 8.5%. The uncertainty 
of the B reco combinatorial background subtraction is esti- 
mated by varying the signal shape function (A r = 3.8%). 
The impact of the binning is studied by changing the bin- 
ning for m x > 1.55GeV/c 2 (A,. = 2.9%). The branch- 
ing fractions of B — > D^*'**Hv and of inclusive and exclu- 
sive D mesons decays are varied within the world aver- 
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Table 1. Fit results for several data samples. 



Sample 


N d 


K 


N c 


«„/../ (*) 


m x < 1.55GeV/e- 


32210 ±233 


167 ±21 


99 ±6 


1.97 ±0.25 


m x < 1.40GeV/c 2 


32210 ±233 


134 ± 19 


64 ±4 


1.77 ±0.25 


m x < 1.70GeV/e 2 


32210 ±233 


191 ±26 


170 ± 11 


2.11 ±0.29 


neutral B reco 


1 1582 ± 133 


76 ± 13 


21 ±3 


2.46 ± 0.43 


charged B reco 


20583 ± 191 


91 ± 16 


77 ±5 


1.68 ±0.30 


Electrons 


18261 ± 173 


99 ± 15 


48 ±4 


2.26 ±0.35 


Muons 


13934 ± 157 


67 ± 14 


47 ±4 


1.66 ±0.36 


9 > 80% 


4491 ±68 


19 ± 7 


13 ± 2 


1.59 ±0.56 


50% < 9 < 80% 


13298 ± 141 


65 ± 13 


46 ±3 


1.85 ±0.37 


9 < 50% 


14122 ± 170 


82 ± 15 


38 ±3 


2.21 ±0.40 



age uncertainties f 1151 (A r = 4.4%). The limited amount 
of simulated events causes an uncertainty A r = 5.1%. 
The uncertainty in the hadronization of the final state of 
B — > X u fv events is determined by measuring R u / S i in bins 
of charged and neutral multiplicities and performing the 
fit using only the nonresonant signal model instead of the 
hybrid model(A r = 3.0%). We also vary the branching 
fractions for charmless semileptonic B decays by one stan- 
dard deviation [ 15|(A r = 2.8%). The fraction of sig- 
nal events with ss contents is varied by 100% for the 
exclusive component and by 30% for the inclusive one [ 
16|(A,. = 3.7%). In the determination of s" el and s u mx we 
allow the nonperturbative parameters to vary according to 
A = 0.480 ± 0.120GeV and X x = -0.300 + 0.105 GeV 2 , 
obtained by scaling the results in I T17I to 0(1 /m?, a s ) in 
order to match the nonresonant MC generator [ 8j. We 
take into account the correlation of -0.8 between A and 
/li(A, = 17.5%). 

We combine the errors related to the detector and the sig- 
nal and background modeling errors quadratically into the 
systematic error and obtain R u / S i = 0.0197 + 0.0027 ± 
0.0023 ± 0.0034, where the errors are statistical, system- 
atic, and theoretical related to the efficiency determination 
and extrapolation to the full mx range respectively. Com- 
bining the ratio R u / t i with the measured inclusive semilep- 
tonic branching fraction of S(B -> Xlv) = (10.87 ± 
O.W(stai) ± 0.300y,s))% [El, we obtain S(B -> XJv) = 
(2.14 + 0.29 ± 0.26 ± 0.37) x 10~ 3 . Using the relation in [ 
15 1 and the average B lifetime of t b - 1.608 + 0.016 ps [ 
1151 we find 

\V ub \ = (4.52 + 0.31 ± 0.27 ± 0.40 + 0.25) x 10~ 3 (5) 

The first error is statistical, the second refers to the exper- 
imental systematic uncertainty, the third gives the theoret- 
ical uncertainty on the extrapolation of R u / S i to the full mx 
range, and the last error combines quadratically the pertur- 
bative and nonperturbative uncertainties in the extraction 
of \V u t,\ from the total decay rate. 

3 Conclusions 

Two different approaches for the extraction of |V„/,| CKM 
matrix element have been presented. The analysis based on 



mx gives currently the most precise determination of \V u b\- 

This is primarily due to specific advantages of this tech- 
nique: large phase-space acceptance and high purity of the 
sample (signal over background ratio ~ 1.7). The two re- 
sults are consistent and they are in agreement with previous 
inclusive measurements i THA . 
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